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■ In the Abehan projection the Yang-Mills theory contains Abehan gauge fields 

■ (diagonal degrees of freedom) and the Abelian matter fields (off-diagonal degrees) 

, described by a complicated action. The matter fields are essential for the breaking of 

the adjoint string. We obtain numerically the effective action of the Abelian gauge 
Q I and the Abelian matter fields in quenched SU{2) QCD and show that the Abelian 

■ matter fields provide an essential contribution to the total action even in the infrared 
region. We also observe the breaking of an Abelian analog of the adjoint string using 
Abelian operators. We show that the adjoint string tension is dominated by the 
Abelian and the monopole contributions similarly to the case of the fundamental 

J> \ particles. We conclude that the adjoint string breaking can successfully be described 

' in the Abelian projection formalism 

o 



o 



X 
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INTRODUCTION 



The mechanism of color confinement in QCD is one of the most important non- 
perturbative problems in the quantum field theory. One of the most promising approaches 
^ ■ to this problem is based on the existence of the dual objects, called monopoles, which are 
condensed in the confinement phase. This approach - known as the dual superconductor 
hypothesis 0] - is realized with the help of the so called Abelian projection Q of SU (N) 



;_i ■ color degrees of freedom to f/(l)^ ^ degrees of freedom. 



The model was shown to be quite successful in explanation of the confinement of the 
fundamental charges such as quarks (see, e.g., reviews j^). Abelian and monopole con- 
tributions to the inter-quark potential are dominant in the long-range region of quenched 
QCD An infrared effective monopole action has been derived in the continuum limit 

after a block-spin transformation of monopole currents 0,0]- It is a quantum perfect action 
described by the monopole currents. The condensation of the monopoles in the confine- 
ment phase was observed in various numerical approaches 0, S|- In the language of the 
monopole currents the condensation irnplies the formation of the percolating cluster studied 
both numerically and analytically jlol |. 

However, this mechanism has a serious problem even in quenched QCD. Although the 
't Hooft scenario describes the confinement of quarks correctly, this scenario predicts also 
the existence of the string tension for the adjoint charges (gluons) in the infrared region. On 
the other hand gluon charges must be screened at large distances due to the presence of the 
gluons in the QCD vacuum. This screening-confinement problem was extensively discussed 
in recent publications 
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The problem of the screening of the adjoint charges in quenched SU{N) QCD has also 
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been discussed in Ref. [12j. The paper provides arguments that the relevant quantity in the 
confinement mechanism is not the Abehan monopoles but the Z{N) center- vortices which 
can explain the screening problem jl^. In our study we pursue a different approach based 
on the dual superconductor model. 

Consider the screening in a confining Abelian model with the charge-two matter fields 
(take, for example, the Abelian Higgs model with the compact gauge fields). The presence of 
doubly charged matter fields screens the confining interaction between the external particles 
with opposite double charges. This happens due to the pair creation from the vacuum at 
certain separations between the external charges. As a result, the potential between the 
particles fiattens at some distances. It should be stressed that the problem is not only to 
explain the fiattening of the potential but also to show the linear behaviour of the potential in 
the intermediate region. On the other hand the charge-one external fields remain unscreened 
in this model. Namely, the potential is linearly rising at large distances. 

The standard model of the dual superconductor in quenched QCD ignores the existence 
of the off-diagonal gluons. However, these gluons have a charge two with respect to the 
Abelian subgroup and they may explain the flattening of the inter-gluon potential which is 
usually studied with the help of the adjoint Wilson loop. On the other hand the introduction 
of the new degrees of freedom - the off-diagonal gluons - should not violate already achieved 
success of the explanation of the quark conflnement in this model. Indeed, quarks have the 
charge one and doubly charged gluons can not screen them^. These and related issues were 
discussed in Ref. l3| for quenched as well as for full SU{N) QCD. 



From the point of view of a realization of the (modified) dual superconductor scenario 
it seems that we have to keep all charge two Abelian Wilson loops in the effective action 
written by the Abelian link fields to reproduce the screening of charge two. Indeed, the 
theory in terms of the Abelian link fields or the Abelian monopole currents alone becomes 
highly non-local if we integrate out all off-diagonal gluon fields after an Abelian projection. 
Needless to say, such an Abelian effective action is useless. The same problem is more serious 
in the real full QCD, since a fundamental charge is also screened in this case. 

The aim of this paper is to calculate numerically the effective action of quenched QCD 
within the Abelian projection formalism. Contrary to previous calculations of this kind we 
include also the doubly charged off-diagonal gluon fields into the effective action and we 
show that their contribution is essential and thus can not be neglected. We also calculate 
correlators of the adjoint Polyakov loops in the Abelian formalism and observe the screening 
of a properly defined potential between static adjoint sources. 

The plan of the paper is the following. In Section |n] we discuss how the screening 
and confinement problem is solved qualitatively in the framework of Abelian dynamics. 
Section IIIII is devoted to the investigation of the Abelian action for the Abelian gauge and 
matter fields obtained by the inverse Monte-Carlo method. In Section IIVI we discuss the 
potential between the adjoint {Q = 2) charges within the Abelian projection formalism. 
We show numerically that properly defined Abelian potential shows screening of the Q = 2 
charges. Moreover, we observe the Abelian and the monopole dominance for the adjoint 
string tension. Our conclusions are presented in the last Section. 



^ However, we may expect a renormalization of the tension of the string spanned between the quarks due 
to the presence of the double charges 
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II. STRING BREAKING IN ABELIAN PROJECTED THEORY 

The partition function of the Abehan effective theory of quenched SU{2) QCD in the 



infrared region may be approximated in the Villain-hke form 



queue 

m 



^^[Jj = j V9 ^ g-j^((rfe+27rn),AZ)(de+27rn))+iQ(e,J) ^ ^-^s^ 

_^ nez:{c2) 
where D is a differential operator 

D^a + ^A-^ + 7A . (2) 

This operator contains local self-interaction term, the Coulomb term described by the in- 
verse Laplacian, A~^, and additional interactions between nearest neighbors. The coupling 
constants d, /? and 7 were calculated numerically in Ref. To simplify notations we use 
the differential form formalism on the lattice jl6|. 

The partition function (PJ) can be rewritten as a string model [l5l |. 

Zq[J] oc Yl e-^(-(^^)-^'^) , (3) 

5a=QJ 

where we have neglected perimeter terms. This model does not contain the dynamical 
matter fields and therefore the string variable a is always closed on the external current J. 
Therefore there is no source for the string breaking in this model. 

Now let us consider the off-diagonal gluons. The Wilson action of quenched SU{2) QCD 

is 

5 = ^^Trt/^,(s), U,,{s) = U,is)U,is + fi)Ulis + u)Utis), (4) 

where U^{s) is the SU{2) gauge field. 

It is convenient to parameterize the SU{2) link variable Ufj,{s) as = c^(s)m^(s) 

where 



COS0,, s tsm0,,(s)e 



ism.(p^\s)e^^^ ' cos<p^(sj 



Q~iSf,{s) 



Here 6, ip and (p are independent variables defined in the regions — vr < 0^{s), V5^(s) < vr and 
< < f. The field 9 behaves as a U{1) gauge field while the field (p corresponds to 

the phase of the off-diagonal gluon field because under an Abelian gauge transformation, 

Q^^^\s) = diag(e^"(^\ e-^"(^)) , (5) 

they behave as follows: 

Of,{s) ^ 9^{s) - df,a{s) = 9f,{s) + a{s) - a{s + ft) , (Pf,{s) ^ (Pf,{s) + 2a{s) . (6) 



The variable (f)^{s) is not affected by the U{1) gauge transformation. After an Abelian 
projection we can integrate this variable out without harming the U{1) content of the model. 
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In order to get an insight of possible forms of interactions between the Abehan gauge and 
Abehan matter fields we replace the averages of cos (f)^{s) and sin by their mean values: 

cos0^(s) (cos0^(s)) = c, sin0^(s) (sin0^(s)) = s, (7) 

where c and s are functions of the SU (2) coupling constant /3. 

As the Abelian projection, we use the Maximal Abelian gauge which is defined by a 
maximization of the functional, 

^ = ^ETr(^3^M(^V3^i(^)) =E(2cosV^(s)-l), (8) 

with respect to the SU{2) gauge transformations, U^^{s) Ufj.{s) = Q{s)U^{s)Q'^ {s + ft). 
The functional (jH)) is invariant under residual U{1) gauge transformations (0). The local 
condition corresponding to maximization (jH)) can be written in the continuum limit as the 
differential equation {d^ + igA^^){Ajj^ — iAj^) = 0. 

The maximization of the functional (jH} corresponds to the minimization of the variable. 
Thus the observation of Refs. made for the mean values 

c ~ 1 , s < c . (9) 

does not come as a surprise. These relations hold in a wide region of the coupling constant /3. 

Following Ref. we rewrite the action of the model (@)) in terms of the variables 9, ip 
and with the help of the definitions (0). Applying Eq. ((7j) to the original action we get 

^TiUf^uis) = cos{Qf,u{s)) - c^s^cos(e^j.(s) - H^^{s) - Cf,u{s)) 

-c\^ cos(e^,(s) + H,^{s) - C^) + c\s^ cos(e^,(s) + H,^{s)) (10) 
+c^s^ cos(e^^(s) - + - C^) + c^s^ cos(e^^(s) - H^) 

+^3" cos(e^,(s) + C^) + cos(e^,(s) - + H,^{s) - 2C^,(s)), 

where we have denoted the U{1) gauge invariant variables as follows: 

e^.(s) = e,{s) + e,{s + fi)- e^{s + z>) - e,{s) , (ii) 

H^,^{s) = 2e^{s) + ip^{s) - ^^{s + fi) , (12) 
Cf,u{s) = ipf,{s) - Lp^{s) . (13) 

The variable 6 is the U{1) plaquette for the gauge field 6, the variable H describes the 
interaction of the matter field ip with the gauge field 9 and the variable C corresponds to 
the self-interaction of the matter field. The validity of the mean field approximation based 
on a self-consistent substitution (|7j) is not known. When we perform the integration, we 
generally get an effective action in terms of Q^u, H^^, and C^^. Below we use numerical 
method to find this effective action. 

A few remarks about the action (llUj) are now in order, (i) From Eq. Q one can immedi- 
ately observe that the leading contribution to the action is provided by the first QED-like 
term depending on the variables 6 only. The coupling between the gauge field 6 and the 
matter field ip is suppressed and the self-interaction of the matter field is suppressed even 
further, (ii) The action (|TU|) should acquire corrections from the Faddeev-Popov determinant 
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resulting from the fixing of the Maximal Abelian gauge. This determinant is an essentially 
non-local functional and the leading local terms were calculated in Ref . jl8l| . 
Let us assume for simplicity the following effective action: 

S,f,. = S^'\e) + S'^'\e,^), = -F,{H)-F2{H')-F,{C). (14) 

where we put H = H^y{s), H' = H^^{s), C = C^iy{s) and Fi, F2, F3 are periodic functions. 
Following Ref. we rewrite the corresponding partition function Z with the external 
source J as follows: 

Scf{.+iQ{9,J) _ / -Dfl-D/o ^-5(i)(e)-S(2)(9,</.)+iQ(e,J) 



Zq[J] = j veVife-^^^-^'^^^^^^ = J VeVipe 



Dip e 



Fi{H)+F2{H')+Fz{C) 



(15) 



The part in the square brackets can be expanded in the Fourier series, 



1=1,2,3 



where n^'^\ i = 1,2,3 are integers and the lattice tensors H, H', n^^\ n^"^^ sum only for n > u 
because H, H' are not anti-symmetric contrary to (C, n^^-*). 

Integrating over (p and summing over n^^-* one can rewrite Eq. (fTBj) as follows: 

[•■■]= Hj)e''^'''\ (17) 

5j=0 

where w{j) are certain weights for the closed current j which is defined from the variables 

n^^^ and n^'^K 

Ms)= E<n^)- (18) 

The general form of Eqs. ()17|) and ()18|) follows from the fact that the fields (f are doubly 
charged and from the gauge invariance of the expression under the exponential function in 
Eq. (fT^. We also give a detailed derivation of Eqs. (fTTj) and (fTSj) in Appendix 1X1 

To simplify further considerations let us rewrite the first term in Eq. ()14j) in the Villain 
form as in Eq. (^. Then we get for the partition function (fT3j) : 

Zq[J] =f'D9j2 Yl ^(•^') ^^p{-^((^^ + 2^^)' ^^(^^ + 2™)) + ^(^' 2j + gj)} . 

5j=0 

Analogously to Eq. Q we get the following model for the string variables dual to the gauge 
field 9: 

Zq[J]= E E Hj)exp{-n'{a,{ADr'a)]. (19) 

Sj=0 Sa=2j+QJ 
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The string model p9|) is different from the model Q by the presence of the doubly 
charged currents representing the contribution of the off-diagonal gluons (the first sum in 
Eq. (|T^). The second sum in this equation is over the integer valued string variable which 
has the dynamical current j as its boundary. 

If the external charge has a unit value, Q = 1, then the dynamical current j can not 
screen the external current QJ and therefore the string always spans on the trajectories of 
the external currents, 6a = 2j + QJ ^ 0. However, if the external current is doubly charged, 
Q = 2, then there exists the dynamical current j = —J such that 6a = 0. This state breaks 
the string: when the distance between the external charges is large enough the state with 
j = — J provides a dominant contribution to the partition function. 



III. EFFECTIVE ACTION FOR GAUGE AND MATTER FIELDS 

In this Section we calculate numerically the effective action for the Abelian gauge and 
the matter fields in quenched SU{2) QCD. We have chosen a trial action in the form: 

S,sie,if) = a,S^ie) + a2S2ie) + a,S,{e) + f]^S,ie,^) , (20) 

where «j, i = 1, 2, 3 and (3i are the coupling constants to be determined numerically. 
The functionals 5*^, z = 1, 2, 3 describe the action of the gauge field 6, 

^1 = - $^[cose^,(s)], (21) 

52 = - ^[cos2e^,(s)], (22) 

53 = +Y1 [si^ ®/^-(^) + /^)] ' (23) 

where the plaquette variable is given in Eq. pij). The action 5*1 is the leading term in the 
Abelian action (fTUI) corresponding to quenched SU (2) QCD in the mean-field approximation. 
The parts S2 3 are also included because they may arise naturally from the integration over 
0. 

As an interaction term between the gauge, 6, and the matter, (p, fields we adopt for 
simplicity 

= -Yl [cos(e^,(s) - H^,{s)) + cos(e^,(s) + H,^{s))] , (24) 

where the plaquette variable H is given in Eq. (fT^ . We have not included other terms from 
Eq. (fTUI) into the trial action because it turns out that the minimal form of the action ()2U|) 
describes the numerical data with a good accuracy. 

We have used the standard Monte-Carlo procedure to generate the gauge field configura- 
tions on the 32^ lattice. The SU (2) coupling constant was chosen in the range /3 = 2.1 ~ 2.7. 
We have generated 100 configurations of the gauge field for each value of the coupling con- 
stant and then used the Simulated Annealing method to fix the Maximal Abelian gauge. 
The couplings ctj, i = 1,2,3 and jSi were determined by solving the Schwinger-Dyson equa- 
tions j20|. We describe the details of this method in Appendix El To make further improve- 
ment of our results towards the continuum limit we used also a blockspin transformation for 
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the SU{2) link variable U^{s): We apply the blockspin transformation to the link variable 

Kis') = ^ I U,{s)U,{s + /,) + 7 ^ U.{s)U,{s + v)U,{s + fi + v)Ul{s + 2/i) | . (25) 

which is visually represented in Figure ^ Here = ^{U) is the normalization factor which 
is introduced to make the fat link belonging to the SU (2) group. The weight parameter 7 
was set to 7 = 0.5. 

FIG. 1: The visualization of the blockspin transformation, Eq. H25|l . 




(c) (d) 

FIG. 2: The parameters a^, i = 1, 2, 3 and f3i for different blocking steps n vs. the scale parameter 
h. The fits by Eq. ()26() are shown by the dashed lines. 

The couplings obtained in this way are depicted in Figures Efa-d). The coupling ai shows 
a perfect scaling since the couphng constant depends only on the physical length 6, (and 
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coupling 


<^o 


C-i 


6o, [tmj 


V 


ai 


0.066(10) 


1.20(2) 


0.61(1) 


2 







0.32(2) 


0.231(7) 


1 


as 





-0.28(3) 


0.46(3) 


1.8(2) 


/3i 


0.064(5) 


0.30(1) 


0.69(2) 


2 



TABLE I: The parameters for the exponential fits (|26j) of the couplings aj, i = 1, 2, 3 and /3i. 



it does not depend on n and a separately). For the couplings 02, «3 and /?i this feature 
does not work: the original data (no blockspin transformation, = 1) is quite different 
from the cases where the blockspin transformation was done (n > 1) while the coupling 
constants with n > 1 scale almost perfectly. One can make a conclusion that the original 
data corresponds to very small values of h where the effective action takes more complicated 
form than (^0)) . 

In order to quantitatively characterize the dependence of the coupling constants on the 
scale factor h we have fitted the data by a function 

/(6) = Co + Ciexp{-(V&o)'^}, (26) 

where Co,i, v and 60 are the fitting parameters. In our fits we have excluded the data without 
the blockspin transformation, ra = 1, for all coupling constants except for a\ case. The best 
fit curves are plotted in Figures El as the dashed lines, and the best fit parameters are shown 
in Table m 

We have found that in the case of a\ and /3i the parameter v is very close to two, and 
therefore in these fits we fixed this parameter, z/ = 2. Similarly, we have also fixed v = \ 
for 02 and Cq = for a2,s- Note that the fit can not describe the coupling a\ accurately at 
small scales, h ^ 0.2 fm. A similar deviation can be found for the coupling as. We expect 
that at small scales the Abelian action becomes much more complicated than the trial action 
(I2(j|21l22l23l24j) which we used to solve the Schwinger-Dyson equations. The small similar 
effect is observed for the effective monopole action obtained by the inverse Monte-Carlo 
methods 1^ • 



The functional 5*1, Eq. (PT|) . makes the leading contribution to the action since the cor- 
responding coupling, tti, is the largest. The actions 5*2 and Ss, in addition to the expected 
action 5*1, play an essential role at small scales since the corresponding couplings, 02 and 
tts, are non-vanishing. The action 5*4, which describes the interaction of the matter fields 
with the gauge fields, has a non-vanishing coupling both at small and large scales similarly 
to 5*1. Moreover, according to Table IHthe couplings a\ and corresponding to these parts 
of the total action, have relatively large lengths compared to the coupling constants 
and as. Thus, at large scales, h\/a 3> 1, the effective Abelian action for the SU{2) gauge 
theory can be approximated as a sum of the QED-like action for the gauge field, Si{6), and 
the interaction term 84(9, ip). 

We interpret the results, obtained in this Section as the manifestation of the Abelian 
dominance (non- vanishing dominant coupling ai) and the importance of the off-diagonal 
(matter) degrees of freedom (non-vanishing coupling /3i). The matter fields are essential 
for the breaking of the adjoint string. From the point of view of further analytical study 
the results of this Section are qualitative because in order to make a quantitative analytical 
predictions at a finite value of the scale b we need much more terms in the trial action fld^ 
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than we have imposed. Indeed, in Ref. [15| the monopole contribution to the string tension 
has been calculated using the effective monopole action. The monopole action was obtained 
numerically and it turns out that in order to get a correct analytical result for the string 
tension one should take into account not only the most local terms in the effective monopole 
action but also a series of the non-local terms. The situation with the effective action for the 
Abelian fields pn|) should be similar to the case of the monopole action since these actions 
are related to each other Nevertheless the adjoint string breaking can quantitatively be 
discussed within the numerical approach on the basis of the Maximal Abelian gauge fixing. 
This topic is discussed in the next Section. 



IV. Q = 2 POTENTIAL FROM POLYAKOV LOOPS 

The easiest way to observe numerically the string breaking effect is to consider the theory 
at finite temperature and define the potential with the help of the Polyakov loop correla- 
tors |ilE3: 



(P(f)pt(^)) = e-v(x-3/)/T^ (27) 



Here T is temperature. 

The adjoint Polyakov loop, Pi, is defined as follows: 




ll^i[U.]j=l{^pl-l), (28) 



where the color vector p = pQ + ip- a defines the fundamental Polyakov loop, P1/2 = l/2Trp, 
p = Hiec^*' ^ straight line parallel to the temperature direction. The adjoint 

Polyakov loop (j2Hl) contains the charged term, Q = 2, and neutral term, Q = 0: 

Pq=2 = liPl- pI) > Pq=o = ^(2pg + 2pl - 1) . (29) 

The Abelian dominance in the most general sense means that a non-Abelian observable 
can be calculated with a good accuracy with the help of the corresponding Abelian operator 
in a suitable Abelian projection. The Abelian dominance was first established for the tension 
of the chromoelectric string spanned between the fundamental sources In this case the 
non-Abelian Wilson (or, Polyakov) loop was replaced by its Abelian counterpart. 



However, in the case of the adjoint potential we immediately encounter a problem [21 
in the Abelian projection the Q = 2 charged component of the Wilson loop shows the 
area law while the neutral Q = component is constant. Therefore, strictly speaking, the 
straightforward Abelian projection of the adjoint operators leads to the vanishing Abelian 
string tension. The simplest way to overcome this difficulty is to introduce the obvious 
prescription for the adjoint operators proposed originally in Ref. 23|. Namely, one should 



disregard the Q = component of the Wilson loop operator and consider the Q = 2 Abelian 
component of the Wilson loop as the Abelian analog of the full (non-Abelian) loop. In 
Ref. 2^1 some numerical arguments in favor of the validity of this prescription were given. 



Below we follow this recipe and show that the string breaking effect can indeed be seen in 
the Q = 2 Abelian and monopole components of the potential. Moreover, we have observed 
the Abelian and monopole dominance for the adjoint string tension. 
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After the Abelian projection the Q = 2 component becomes 

P^t2 = cos2^c, ^c = $^^.. (30) 



where 'dc enters the Q = 1 Abehan Polyakov loop, P^^^ — cos'dc- 

We calculate numerically the static potential between the adjoint particles using the 
Polyakov loop correlators (jTTj) . We use four types of the Polyakov loops: non- Abelian, 
Abelian, monopole and photon Polyakov loops, 

Pq=2=pI-pI Pa=2 = cos2^c, P<^°^ = cos2^--, P^l, = cos2i}f , (31) 
respectively. 

The functions 19™°° and represent the contributions to the Polyakov loop coming from 
the monopole currents and the photon fields, respectively 0,13: 

= - 5^ 5^ I^(f - f, t - t')dle,,{x\ t') , (32) 

t x',t' 

= _27r J] ^Z}(x-x',t -0-9X4(5', 0, (33) 



t x\t' 



where the variables O G (— vr, tt) and n E Z are extracted from the Abelian plaquette 
variable, Q^y{s) = 6*^(5) + 6^{s + fi) - d^{s + 0) - 6*^(5) = Q^y{s) + 27m^^{s). D{s) is the 
inverse Laplacian, d'^d^D{s) = —Sq^s- 

We numerically measured the potential between the static adjoint sources on the 16^ x 4 
lattice at /5 = 2.2 (confinement phase) using 2000 configurations. The Abelian, monopole 
and the photon components of the potential were measured in the Maximal Abelian gauge. 
In order to reduce the statistical errors in our calculations of the potentials we have ap- 
plied the Hypercubic Blocking 2^ procedure to ensembles of the non-Abelian, Abelian 



and photon gauge fields. We have not applied the blocking to the monopole contribution 
of the potential because in this particular case the blocking makes the data noisier. The 
Hypercubic Blocking method is briefly described in Appendix O 

We present the numerical results in Figure |21 One can clearly see that all potentials 
become flat in the infrared region clearly indicating the presence of the string breaking. 
The non-Abelian potential as well as the Abelian and the monopole contributions contain 
the linear pieces at small enough distances while the photon contribution to the potential 
does not contain a linear part. These observations are in a qualitative agreement with the 
Abelian (monopole) dominance hypothesis 

To make a quantitative characterization of the potentials we fit our data by a function 

exp| — |=exp| |+exp| 1, (34) 

where we have chosen the string potential in the simplest form, V^tr(-R) = cfq=2R- The fitting 
parameters are the adjoint string tension o"q=2, the mass parameter m and the self-energy 
Vq. The first term in Eq. (j34p corresponds to the broken string state and the parameter m 
is the mass of a state of "external heavy adjoint source" -"light gluon". The second term is 
the unbroken string state. Here we neglect other states including the string excitations. 
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FIG. 3: The potential between adjoint static sources and the Abelian, the monopole and the photon 
contributions to it. The fits by the function (|34() are shown by the solid lines. 



type 




m/y/a 


Non-Abelian 
Abelian 
Monopole 


2.49(3) 
2.33(3) 
1.94(1) 


1.28(1) 
1.84(1) 
2.27(2) 



TABLE II: The parameters for the fits of the potential by the function (|34)) . Here a = aii2{T = Q). 



We perform the fits in the range starting from two lattice spacings, rmin = 2a. The reason 
for this restriction is twofold: (i) the hypercubic blocking modifies the potential at small 
distances; (ii) in our fitting function (jMj) the perturbative Coulomb interaction (which is 
essential at small distances) is not included^. 

The best fit functions are shown in Figure El by the dashed lines and the best fit param- 
eters are presented in Table lUl One can clearly see the existence of the Abelian dominance 
for the string tension: cXgLg ~ 0.94 (Tq=2, where (Jq=2 is the string tension extracted from 
the non-Abelian Polyakov loop correlator. The monopole dominance can also be observed: 
^Q=2 ~ O.SSag'^g ~ 0.78 (Tq=2. The monopole dominance is less manifest than the Abelian 
dominance in agreement with the precise observations at /? = 2.5115 in the case of funda- 
mental external sources, Ref. 

In Ref. jsf the potential between the static Q = 2 Abelian sources has been measured in 
the zero temperature case. Despite the string breaking has not been observed in this case. 



^ Nevertheless we have checked the effect of the Coulomb interaction shifting the string potential as 
14tr(-R) Veti-iR) — a/R, where a is an additional fitting parameter. We have observed that the best 
fit values of the parameters (Tq=2 and m got the shift about 1-2% which is of the order of the statistical 
errors for these parameters. 
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the ratio between Q = 2 and Q = 1 Abelian string has been measured: (Tq=2/c"q=i = 2.23(5). 
Taking into account that the ratio between Q = 1 Abehan and SU{2) string tensions is jsf, 
(Tq=2/o- = 0.92(4), we get the prediction of Ref. jl] for the ratio 0-^=2/0" = 2.42(12). We 
observe a very good agreement with our resuh, aQ=2/<^ = 2.33(3), given in Table ITTl 

According to our numerical results the Abelian and the monopole contributions to the 
masses of the heavy-light adjoint particles, m, do not coincide with the corresponding mass 
measured with the help of the non-Abelian Polyakov loops. On the other hand we do not 
expect neither Abelian nor monopole dominance to hold in this case since these types of 
dominance are usually valid for infrared (non-local) quantities in accordance with the ideas 
of Ref. 0. Because of the local nature of the mass m the Abelian/monopole dominance 
may not work in this case. 

The absence of the Abelian dominance for the mass parameter m implies the absence 
of the Abelian dominance for the string breaking distance. Indeed, the simplest definition 
of the string breaking distance, Rsh, corresponds to a value of R at which both terms in 
Eq. are equal. For the hnear string potential, Vstr = o'q=2R, this distance is defined 
as Rsh = 2m/(jQ=2- In other words, the string breaking distance is the distance where the 
energy of the string, o"Q=2-Rsb, is equivalent to the energy of the two heavy-light states, 2m. 
Since the Abelian dominance works only for the string tension a"Q=2, the string breaking 
distance, Rsh, should not be Abelian/monopole dominated quantity. 



We have calculated the effective action for the Abelian gauge and the Abelian charged 
matter fields in the Maximal Abelian projection of quenched SU{2) QCD. We have shown 
that in the infrared limit the contribution of the matter field to the action is non-vanishing. 
Thus we have shown at the qualitative level that the matter fields, carrying the Abelian 
charge Q = 2, must lead to the adjoint string breaking. To check this effect on the quantita- 
tive level we have studied the potential between adjoint static sources as well as the Abelian 
and the monopole contribution to this potential. We have observed the string breaking 
(flattening of the adjoint potential) manifests itself in the Abelian and the monopole contri- 
butions similarly to the non-Abelian case. Moreover, we show that the adjoint string tension 
is dominated by the Abelian and the monopole contributions analogously to the case of fun- 
damental particles. Thus we conclude that the adjoint string breaking can qualitatively be 
described in the Abelian projection formalism. The key role in the adjoint string breaking 
in the Abelian picture is played by the off-diagonal gluons which become the doubly charged 
Abelian vector fields in the Abelian projection. 



V. CONCLUSIONS 



APPENDIX A: DERIVATION OF EQ. (fT7|l 



In this Appendix we present a detailed derivation of Eq. (fTTjl . 



n 




Sj=0 



(Al) 



The Fourier transformation applied to each of the terms in the l.h.s. of this equation 

"3 



Vif e 



Fi{H)+F2{H')+Fz{C) 



1=1,2,3 



where /j are the Fourier components of e^' and $ is the phase: 



.IJ,>U 



Using definitions (fTT| - (P|) we get: 

The integration over the field if gives two constraints 



77 

which lead to the 



<^ is) - - A) + <J is) - ^ is-i>)=0. 

Equation ()A4|) gives the natural definition for the current of the matter fields: 



j,is)= E<^(^)- 

v(<At) !^(>M) 

Note that due to the constraint (IA6I) the current (IA7I) is closed: 



Combining Eqs. (jX2| . (|X4|) and (jXTj) . we get the r./i.s of Eq. (jXTJ with 



j=1.2,3 



s i=l 



i/(</i) 



s — u 



+ 2n!f)(.; 
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APPENDIX B: SCHWINGER-DYSON EQUATIONS 



Consider a model of the gauge field 9. The expectation value of an arbitrary operator 
0{6) measured at the ensemble {6i} of the gauge fields 6 is 



{0{e)) = j veo{e)e~''-'^ = l[j de,o{{e,}) e 

Shifting one of the hnk fields Qi^ at the link zq by an infinitesimal value e we get 

TT TT TT 



(Bl) 



TT TT + e 



TT TT 



'9^ 



fB2) 



The requirement that this shift does not change the partition function gives the Schwinger- 
Dyson equation: 



7r TT 



n /rf^J {0{K.m.^^.) e-^(^-^^^>...)] = / [0(^)e"^(^)] = 0, (B3) 



which can also be rewritten in the form: 

/dO{e)\ 



ds{e)\ 







(B4) 



To determine the parameters of the trial action (j!2(J|) - (j24j) we solve Eq. ()B4|) with the 
following set of operators: 



Oi 



dSi 



(1 = 1,2,3,4), O5 



(B5) 



The expectation values of these operators give a set of five Schwinger-Dyson equations: 

3 



5^ a J 

.7=1 



I3i 



dSj dSj 

2^ 



Pi 



dSi dS, 



J 



(1 = 1,2,3,4), (B6) 
(B7) 
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Since we have five equations ()B6|) - ()B7|) to determine four independent couplings, ctj, 
i = 1,2,3 and Pi, the system of equations ()B6|) - ()B7|) is overdefined. Thus we find the 
couplings with the help of Eq. ()B6|) and then use Eq. ()B7|) as a consistency check. We 
find that for the original fields the l.h.s. of Eq. ()B7|) is approximately 10% larger than the 
r.h.s. However, after applying the blockspin transformation the discrepancy becomes much 
smaller (it becomes of the order of the statistical errors), and the solution of Eqs. ()B6pB7jl 
becomes self-consistent. 



APPENDIX C: HYPERCUBIC BLOCKING 



The Hypercubic Blocking (HYP) procedure is a version of the smearing method which 
allows to reduce the noises of the lattice gauge fields jl^ . As a result the statistical errors 
of ensemble averages of various operators are reduced. HYP is replacing gauge link fields, 
U^{s), by "fat links", V^(s), according to the following scheme: 



1 

h 
1 

h 
1 

h 



ai)Uf,{s) 



ai 



(CI) 



as 



where ki, i = 1,2,3 are chosen in such a way that the matrices ()C1|) belong to the SU(2) 
group. We choose the parameters of the HYP, ai, 0:2, 03 G [0,1], following Ref. j23[: 
«! = 0.75, a2 = 0.60 and 03 = 0.30. At these values the smoothing of the gauge field 
configurations is most efficient. 
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